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Spins transverse to the magnetization of a ferromagnet only survives over a short distance. We 
develop a drift-diffusion approach that captures the main features of transverse spin effects in systems 
with arbitrary spin textures (vortices, domain walls) and generalizes the Valet-Fert theory. In 
addition to the standard characteristic lengths (mean free path for majority and manority electrons, 
spin diffusion length) , the theory introduces two lengths scales, the transverse spin coherence length 
£± and the (Larmor) spin precession length £l- We show how £l and £± can be extracted from 
ab-initio calculations or measured with giant magneto-resistance experiments. In long (adiabatic) 
domain walls, we provide an analytic formula that expresses the so called " non-adiabatic" (or field 
like) torque in term of those lengths scales. However, this "non adiabatic" torque is no longer a 
simple material parameter and depends on the actual spin texture: in thin (< lOnm) domain walls, 
we observe very significant deviations from the adiabatic limit. 

PACS numbers: 



Soon after the seminal paper of J. Slonczewski which 
computed the spin transfer torque [T] in a non collinear 
spin valve, the concept of transverse spin current be- 
came widely discussed in the spintronics community. The 
problematic of transverse spin current can be formulated 
in a simple setup: using a first ferromagnet, one injects a 
spin polarized current into a second ferromagnet whose 
magnetization is perpendicular to the first. A natural 
question that arises is what happens to this spin polar- 
ization which is transverse to the magnetization? One 
can simply surmise that after a short distance the polar- 
ization of the current becomes aligned with the magneti- 
zation of the second magnet. The corresponding loss of 
transverse spin current is assumed to be transferred to 
the magnetic degrees of freedom, hence a spin torque is 
exerted on the magnetization. In spin valves, a good un- 
derstanding of the physics could indeed be reached by as- 
suming a full absorption of the transverse spin current at 
the interface, i.e. a purely interfacial spin torque [5] ■ In- 
deed, quantum calculation soon showed [3] that the lead- 
ing mechanism for the absorption of transverse spin cur- 
rent is of purely ballistic origin: the band structure be- 
ing in general very different for minority and majority 
spin in a ferromagnet, the resulting fast spatial preces- 
sion of the transverse spin leads, upon averaging on the 
incident directions, to an exponential decays of the trans- 
verse spin current as a function of the distance of pen- 
etration in the ferromagnetic layer. The transverse spin 
coherence length £± obtained from these calculations is 
typically rather small, a few nm, which justifies the inter- 
facial limit. There exist, however, many situations where 
the interfacial limit £j_ — > entirely misses the relevant 
physics. One of those cases is current induced domain 
wall motion. It was recognized early that, in domain wall, 
the main " interfacial" spin torque (known in this context 
as the "adiabatic torque") is not sufficient to understand 



how a current can set a domain wall in motion[H [5]. A 
second spin torque, perpendicular to the main one and 
usually much weaker (known as the "non adiabatic" of 
"field like" torque) is necessary to describe the dynamics. 
An important theoretical effort has been devoted to the 
calculation of this "non adiabatic" torque torque using 
a wide variety of techniques ranging from quantum [BJ- 
12] to phenomenological {10 approaches. This effort also 
include a large variety of experimental works |llH15j 

In this letter, we perform four things. First, we de- 
velop a drift-diffusion theory capturing finite transverse 
spin current effects in systems involving ferromagnetic 
diffusive metal regions with arbitrary three dimensional 
texture. The theory can be thought as a direct gener- 
alization of the Valet-Fert theory[16] to arbitrary non 
collinear systems beyond a two current formulation in 
the ferromagnetic regions. In the lumped circuit element 
(discrete) limit, it is equivalent to the generalized cir- 
cuit theory [TTJ [TB]. Second, we show how the two new 
parameters of the theory (spin coherence length £± and 
the (Larmor) spin precession length £^) can be obtained 
from ab-initio quantum calculations. Third, we propose 
an experimental setup that allows to measure l\_ and £l 
directly using current perpendicular to plane (CPP) gi- 
ant magneto-resistance (GMR) measurements. Fourth, 
we apply our theory to spin torque in domain wall hence 
providing a direct link between CPP GMR and the "non 
adiabatic" torque. 

Generalized Drift Diffusion Theory. Our starting point 
is a set of equations for the current densities in the charge 
j^{r) and spin j a (r) sectors and local charge fJ, c (r) and 
spin \x chemical potentials of the system. The entire the- 
ory is fully equivalent to the Continuous Random Ma- 
trix Theory (CRMT)pIl EE] developed by some of us 
and is obtained from the latter through a simple change 
of variable after some straightforward, though somewhat 
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lengthy, calculations 29 . The actual details of the deriva- 
tion will be postponed to a subsequent publication and 
we concentrate here on the physics implications. Simi- 
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where the explicit index a G {x, y, z} stands for the spa- 
tial direction with d a = d/da while bold vectors are used 
for the three dimensional spin space. The unit vector 
m(r) points along the local direction of the magnetiza- 
tion. The "charge" and "spin" currents defined above 
have the dimensions of energy. They are simply related 
to the (observable) electrical current density, 

I a = 4fJ(en Sh ) (5) 

and spin current density 

J a = 2hjJ(e 2 TZ sh ) (6) 

where IZsh is the Sharvin resistance for a unit surface 
(typically 0.5/f2.m 2 ) and e < the charge of the electron. 
The theory is parametrized by five independent lengths 
scales: the mean free path for majority £^ and minority 
£^ electrons, the spin flip diffusion length £ s f, the spin 
coherence length £± and the (Larmor) spin precession 
length £l. Alternatively, one can introduce the average 
mean free path (l/£* = + l/^j.) and polarization 
ft = (if — l±)/(l\ + These parameters are totally 
equivalent to the usual parameters of Valet-Fert theory 
[iniHH] with the following correspondance: i a — TZgh/pa 
(p a is the spin dependent resistivity) , same (3 and £ s f and 
4 = ft sh /(4p*)- 

The physical meaning of £± and £l is best identified 
by studying the transverse spin in a non textured magnet 
d a m = 0. n = /Lin + H± (and j a ) is decomposed into 
a longitudinal fi^ = (/x • m)m and perpendicular /x^ = 
(m x fx) x m contribution. Equations Q-Q reduce to 
Valet-Fert equations for the charge and longitudinal spin 
part. Expanding fi = fj,i\m + [i\e\ + ^2^2 with e\ and 
e2 two orthonormal unit vectors perpendicular to m we 
arrive at J2 a daafi = AAmx for M = A*X + ^2 with 




larly to Valet Fert, the set of equations consists of gen- 
eralized Ohm laws and current conservation equations: 



(1) 
(2) 
(3) 

(4) 
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which for most systems (except Nickel, see below) reduces 
to l/l mx ~ — i/t-L- In other words, the transverse 
spin accumulation jl(x) cx exp(— x/l mx ) decays over a 
length scale £± and precesses around m over a length 
£l. We note that other authors have proposed general- 
ization of the drift-diffusion equations before [TU1 |2"UI |2"T] . 
While these approaches captured the precession part of 
the above theory, they suffer from the absence of the 
(crucial) terms with £± so that the role of absorbing the 
transverse spin is taken by the (much larger) lengths £„ 
and £ sf (Eq.Q with £±=oo). 

The equations for a non magnetic metal are obtained 
from Equations ([lJ-Q by setting £j_ — > 00 and (3 = 0. 
The presence of a magnetic field B (at the origin of the 
Hanle effect) is captured using £l = hvp / (g^BB) and 
m = B/\B\. 

Interface and reservoirs boundary conditions. To com- 
plete the theory, we need the boundary conditions be- 
tween one normal (n) and one Ferromagnetic (f) mate- 
rial. Noting n a the vector normal to the interface and 
pointing toward the magnetic material, A// = \x n — fj,^ 
the difference of chemical potential across the interface, 
e n = — ef = 1 and a = n/ /, we get: 

= °* i( m ■ A M) + 7A/i c ] m + 3f?« nx )A/Xj_ 
- 3« x )m x A/x_l + e a [9?« X K - 3« x )m x 

(8) 

n af a = G * l A Vc + jrn, ■ Afi] (9) 



where 7 is the (Valet-Fert) polarization of the inter- 
face resistance and a* is related to the Valet-Fert as 
1/rr* = 2rjJ(l — 7 2 ). The other "mixing" parameters are 
expressed in term of the mixing transmission (T mx ) and 
reflection (i? mx ) parameters [IB] of the interface as fol- 
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lows, 
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Last, the boundary conditions between a normal elec- 
trode at a potential eV and the system reads (n Q points 
toward the system), 



Ua ^ + M = 

a 



(13) 
(14) 



Extracting l\_ and £l from ab-initio calculations. 
While an important experimental work has been devoted 
to the calibration of /3, and £ s f 7 very little is known 
about the actual values of the lengths associated to trans- 
verse spins. A first insight is given by ab-initio calcu- 
lations which measure the mixing transmission T mx of 
simple Normal- Ferromagnet-Normal metal trilayers |23j . 
T mx is a complex number whose amplitude measures the 
probability for a transverse spin to go through the system 
while its phase measures the angle of precession. Up to 
tiny corrections (due to multiple reflections of the trans- 
verse spin at the interfaces), it is given by, 



(15) 



where d is the thickness of the magnetic layer and T™* 
the mixing transmission of the Normal-Ferro interface. 

and 



Eq.(15l allows for a direct extraction of £± and £^. An 
example of ab-initio data, taken from Ref. [23 is shown in 
the upper part of Fig [l] where we plot the real and imagi- 
nary part of T mx for a Cu-Co-Cu trilayer (the inset shows 
the phase). We find that a good fit with Eq. (15) could 
be obtained allowing to extract £± and £l with good 
precision (the value of the interface mixing transmission 
being somewhat less accurate). We have repeated the 
procedure with ab-initio data available in the litterature 
and collected the results in Table |TJ We find typical value 
lh ~ 0.3nm (corresponding to a full precession of 2-7T on 
2nm) and £± 2nm. Nickel, a somewhat weaker magnet 
seems to have a significantly longer transverse coherence 
length £± ss 6nm than the other materials. We note that 
since £j_ and £l originate from ballistic effects, they de- 
pends, in principle, of the crystalline direction as seen 
in Table [I] Such an effect could be incorporated into our 
theory by using tensorial instead of scalar values for these 
lengths. However, given the lack of information on these 
lengths, we restrict to scalar values at this stage. 

Measuring £l and £± from CPP GMR. An alternative 
route to the ab-initio calculations is to actually measure 




FIG. 1: Upper panel. Cu-Co-Cu trilayer. Real (circles) 
and imaginary (diamonds) parts of T mx obtained from the ab- 
initio data of Ref. 22 as a function of the thickness d of t he Co 
layer. The lines are the corresponding fits with Eq. ( 15 1. Left 
inset: schematic of the setup. Right inset: same as the main 
panel for the phase of T mx . Lower panel. Inset: schematic 
of the proposed experimental setup for the measurement of 
£± and £l- the studied magnetic layer X is sandwiched be- 
tween an analysing A and polarizing layer P with normal 
spacers S. The studied layer must have its magnetization 
perpendicular to the one of A and P. Main plot: Numerical 
calculation of the GMR (in % ) as a function of the thick- 
ness d of the studied layer for a stack with S = Cu and 
A = P = (Cuoa | Nio.s)x3 (multilayer with perpendicular 
anisotropy). The squares (circles) show the numerical data 
for X = iVi(lOO) (X = Co(lll)) while the lines correspond 
to the fit with Eq.fili]). 



£l and £±, in the same spirit as what has been done for 
the other three "Valet-Fert" parameters in the collinear 
configuration. We propose the setup shown in the lower 
panel of Fig[T] which consists of a standard spin valve in 
the middle of which one inserts the layer that one wants 
to study (X). The magnetization of the X layer must be 
perpendicular to the ones of the spin valve (A and P) so 
that the latter must have strong perpendicular anisotropy 
if X has planar anisotropy or vice versa. One performs 
standard GMR measurement and measures the resistance 
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in the parallel (Rp) and anti-parallel (R ap ) configuration. 
The GMR = (R ap — Rp) / '{R ap + R P ) signal is proportional 
to the real part of the mixing transmission of X (and its 
interface) so that the expected signal reads, 



GMR(d) = A cos 



d 



(16) 



where the two constants A and S depend on the material 
parameters of the spin valve. An example of the expected 
signal is shown in Fig [I] where we chose (Cu .4 | Nio.s)x3 
(indices are in nm) as our polarizing and analyzing layers 
with perpendicular anisotropy |26j . The symbols show 
the numerical calculations for X = Ni and X = Co 



together with the fit with Eq.( 16 ). We find that the GMR 



signal, though a bit smaller than in standard spin valve, 
lies around 1% with a very clear oscillating pattern. 

Definition of spin torque. Before turning to a practical 
calculation of spin torque in a domain wall, we need to 
identify its proper definition. In the original work of Slon- 
czewski [T], spin torque was defined using a very robust 
conservation argument: whatever spin current has been 
lost by the conducting electrons must have been gained 
by the magnetic degree of freedom (conservation of to- 
tal spin). This is not true in the presence of spin-orbit 
coupling however (which is chiefly responsible for the fi- 
nite £ s f in metals), as part of the angular momentum 
is transferred to the lattice. Hence, in the spin conser- 
vation equation Q, the divergence of the spin current 
is the sum of a (spin-orbit induced) spin flip term and 
the term corresponding to the exchange coupling to the 
magnetization. The latter corresponds to the spin torque 
density t which reads, 



2h 



e 2 K 



Sh 



jj- (™ x fi) x m - j- (m x n) ) U 7 ; 



Field like torque in domain walls. We now turn to a 
system with a non trivial magnetic texture and study spin 
torque in a one dimensional domain wall. We redefine 
the local basis as &\ — d x m/\d x m\ and e 2 = m x e x . 
Noting 9(x) the angle between m(x) and the z axis (with 
9 = d x 9) we obtain, 



dxxfi + d x (n\\6) = 
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We now proceed with the " adiabatic" limit and consider 
the case of a very long domain wall (so that, up to small 
corrections, the spin accumulation follows adiabatically 
the magnetization[2"T] ) . We obtain, 
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(19) 



from which we can calculate the torque. The torque T\ 
along ei is the main contribution predicted by Berger[2"S] 



(each up electron leaves the system as a down elec- 
tron and hence deposit h on the domain wall) T\ ~ 
—f39HI/(2e). The torque t 2 along e 2 is the "non adia- 
batic torque" (or "field like" torque) whose importance 
has been stressed in the introduction. Introducing the so 
called "beta" parameter, f3 T = — t 2 /ti, we arrive at 



IjJU 
P 2 



(20) 



The importance of Eq.(20) comes from the fact that it 
connects two different realms: domain wall motion on the 
left hand sight and CPP GMR physics on the right hand 
side. In the limit of long spin-flip lengths, Eq.(20) re- 
duces to /3 T rj £l^*I '^ s f (up to a prefactor of order unity) 
which is similar (but not equivalent) to the one obtained 
by Zhang and Li[TU] (3 T as i^/igf. However, in the limit 
of strong spin-flip scattering, the parameter j3 T saturates 
toward j3 T ~ £±/£l. The different values of this param- 
eters have been summarized in Table [I] Equation. ( 20 ) 
seems to imply that f3 T is a simple combinations of ma- 
terial parameters, but it is only valid for very smooth 
variation of the magnetization. Fig(2] studies numerically 
the parameter f3 T at the center of the domain wall as 
a function of the width lyy of the wall (a simple form 
tan[0(x)/2] = exp(x/lw) was assumed for the wall). We 
find that the adiabatic limit Eq. ( 20 ) works very well for 



wall thicker than lOnm but a very significant increased 
is observed for thinner walls. Values close to unity are 
expected for very thin walls (such as those found in sys- 
tems with strong perpendicular anisotropy) and strong 
spin-orbit coupling. 

Conclusion. The drift-diffusion approach developed 
here has the important advantage of hiding many mi- 
croscopic details which eventually only leads to a renor- 
malization of the parameters of the theory. On the other 
hand, it does capture the crucial physical ingredient of 
transverse spin physics: its rapid absorption due to de- 
coherence between different directions of propagation. 
Once the effective parameters of the theory are measured 
(as it was done in the collinear configuration) or calcu- 
lated, a large number of predictions and links can be 
made. For instance, CPP GMR and domain wall mo- 
tions are usually considered as involving quite different 
physics, but here we have shown that a direct connexion 
can be made between both. 
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TABLE I: Material parameters for various materials. The reference corresponds to the ab-initio data from which £l, ix and 
TJ-Ux have been extracted. A unique value of the Sharvin resistance of TZsh = 0.5/fi.m 2 has been used everywhere. The domain 
wall "non adiabatic" j3 T parameter has been computed using Eq. (201. 



